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r^, Abstract 

^""i We investigate the possibilities of New Physics affecting the Standard Model (SM) 

Qs Higgs sector. An effective Lagrangian with dimension-six operators is used to capture the 

^O effect of New Physics. We carry out a global Bayesian inference analysis, considering the 

*K^ recent LHC data set including all available correlations, as well as results from Tevatron. 

—J Trilinear gauge boson couplings and electroweak precision observables are also taken into 

^^ account. The case of weak bosons tensorial couplings is closely examined and NLO QCD 

corrections are taken into account in the deviations we predict. We consider two scenarios, 
one where the coefficients of all the dimension-six operators are essentially unconstrained, 
and one where a certain subset is loop suppressed. In both scenarios, we find that large 
deviations from some of the SM Higgs couplings can still be present, assuming New Physics 
arising at 3 TeV. In particular, we find that a significantly reduced coupling of the Higgs to 
the top quark is possible and slightly favored by searches on Higgs production in association 
with top quark pairs. The total width of the Higgs boson is only weakly constrained and 
can vary between 0.7 and 2.7 times the Standard Model value within 95% Bayesian credible 
interval (BCI). We also observe sizeable effects induced by New Physics contributions to 
tensorial couplings. In particular, the Higgs boson decay width into Z7 can be enhanced 
by up to a factor 12 within 95% BCI. 
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1 Introduction 

The existence of a scalar resonance compatible with the SM Higgs boson has now been 
firmly established at the LHC by the ATLAS and CMS experiments [1, 2]. Based on the 
latest results presented at the Rencontres de Moriond 2013, the combined significance of 
the excess around 125 GeV reaches more than 7 a in both experiments. The analyses are 
based on up to 5 fb'^ at 7 TeV and 21 fb"^ at 8 TeV, collected in 2011 and 2012. 

Following this fundamental breakthrough, new questions need to be addressed. Some 
of them may find an answer during the LHC era, through the analysis of the properties 
of the new boson. A first question is whether or not this resonance is a Higgs boson, i.e. 
the manifestation of a field involved in the electroweak symmetry breaking (EWSB) and 
unitarization of WW scattering. Second, if it is indeed a Higgs, one may wonder to what 
extent it is compatible with the Standard Model Higgs. Indeed, the properties of the Higgs 
boson could be dramatically modified with the theoretical prejudice that New Physics has 
to emerge near the electroweak scale. On the other hand, direct searches for new states 
beyond the SM have so far turned out to be unsuccessful, and indirect constraints from 
electroweak precision measurements at LEP push the limits on masses of new particles 
somewhat further above the electroweak scale. If New Physics is indeed present and is 
somehow separated from the electroweak scale, the couplings of the Higgs boson will be 
close to those of the SM and will only be modified by the effect of a few higher dimensional 
operators. In this paper, we will explore an effective field theory (EFT) with only relatively 
few new parameters. 

Many aspects of these fundamental questions have already been investigated in several 
works (see e.g. [3-12] for studies based on the 8 TeV results). For example, the chiral 
EW Lagrangian with a non-linear realization of the SU{2)l x U{1)y symmetry is a quite 
general framework in order to study the properties of electroweak symmetry breaking. 
Many scenarios producing possibly large deviations with respect to the SM Higgs properties 
are successfully captured in such EFT approach. 

We are going to consider that new states appear at a typical scale A substantially larger 
than the electroweak scale. For physical processes involving an energy scale smaller than 
A, New Physics can be integrated out. As a consequence of this hypothesis, the resulting 
low-energy effective theory consists in the Standard Model, supplemented by infinite series 
of local operators with higher dimension, which involves negative powers of A, 

/:efr = /:sM + E S^- (1-1) 

i 

The effects of such higher dimensional operators (HDOs) have been investigated in many 
contexts such as flavor physics, or the study of the properties of the electroweak gauge 
bosons through LEP precision measurements. The purpose of this paper is to study the 
electroweak sector again, which now includes new Higgs observables. For our analysis, 
we only have to consider the leading HDOs. The only operator with n^ = 1 is the one 
giving a Majorana mass to the neutrino, and is not relevant for our study. We will thus be 
exclusively interested in the n = 2 terms, i.e. dimension-6 operators. 



We adopt the statistical framework of Bayesian inference in this work, which allows us 
to assign probabilies to our parameters and to deal with partially constrained problems. 
Another interesting property is that the unnatural {i.e. fine-tuned) character of precise 
cancellations which may occur between HDO contributions is built-in in this framework. 
Indeed, regions of the parameter space in which precise cancellations occur have by con- 
struction a weak statistical weight (see also [14]). The results we will present can thus be 
considered as generic, i.e. free of improbable cancellations. Our Bayesian analysis will rely 
on Markov chain Monte Carlo (MCMC) techniques, which allows us to easily sample the 
posterior probability distribution function. 

The outline of the paper is as follows. In Section 2, we lay out the formalism for higher- 
dimension operators in the electroweak sector. In Section 3, we present the dataset used for 
the analysis and the measurements entering the likelihood functions. The peculiar case of 
observables sensitive to tensorial couplings relating Higgs and weak bosons is investigated 
in Section 4. In Section 5, we derive the observable deviations from the SM induced by 
the higher dimensional operators, taking into account leading NLO QCD effects. Section 
6 presents the setup of our Bayesian analysis. Section 7 is devoted to our results. Our 
conclusion is given in Section 8. 

2 Electroweak higher-dimension operators 

In this section, we define the basis of dimension-6 operators supplementing the renor- 
malizable electroweak sector of the SM Lagrangian. We refer to [15, 16] for further details 
on the Standard Model HDOs. 

A basis of CP-even operators not involving fermions can be chosen as ^ 



Any other operator can be reduced to these via integration by parts and the use of the SM 
equations of motion for the Higgs and gauge fields, possibly generating operators involving 
fermions. Amongst the latter, only a limited set will be relevant for our purpose. Operators 
of the form Jh ■ Jj, where Jh and Jf are SU{2) or U{1)y currents involving Higgs field 
and fermion / respectively, will in general contribute to FCNC as well as electroweak non- 
oblique corrections {e.g., non-universal couplings of fermions to gauge bosons). However, 
the operators 

OD = J^^J^ Od = 4m^J, (2.4) 

where J = ^f Jf are the SM fermion currents coupling to B^ and VF^, are flavor diagonal 
and only result in universal corrections to gauge couplings and should hence be viewed as 



Hf, 


On2 = \H\^\D^H\\ 


0'^, = \H^D^H\\ 


(2.1) 


ra \2 


Obb = H^H {B^,f , 


OwB = H^W^,HB,,,, 


(2.2) 




OGG = H^H{G%f. 




(2.3) 



^The operator Oe plays no role in what follows and is listed here only for completeness. 



contributing to S and T. ^ We will also need to consider Yukawa corrections of the form 

Of = 2yf\H\^HfLfR, (2.5) 

where fjt = ta, bji, tr and fi the corresponding doublet (/£ = e'^^q^ , g^' i ^l ) ^'^^ ^Z 
the Yukawa coupling. 

Note that the operators Od and O'j-, could be traded for the operators 

Ow = {D^H)^W^,D,H , Ob = {D^H)^D,H B^, , (2.6) 

by use of the SM equations of motion for B and W. While Od and O'j^ contribute to S 
and T but not to the modified Higgs couplings, for Ob and Ow it is the other way around. 
Both choices of basis are physically equivalent. Before passing from a general redundant set 
of operators to a convenient irreducible basis via the equations of motion, it is useful to first 
identify the operators that cannot be generated at tree-level. This is valuable information 
and we would like to avoid it to be lost in the course of the reduction. However this is what 
would happen if we eliminated Oo and O'^ in favor of Ow and Ob- Indeed, this would 
cause the coefficient of e.g. OwB (which cannot be generated at tree-level) to be shifted 
by the coefficient of Od (which can be generated at the tree-level via exchange of spin-one 
states). This is why we choose this basis. 

The only remaining two-fermion operators are of the dipole type. These operators 
are tightly constrained by FCNC as well as by their contributions to electric and magnetic 
dipole moments. Moreover, they are necessarily generated at the loop-level, and only affect 
Higgs couplings to gauge bosons by modifying existing SM loops. They will not have any 
impact on our results, therefore we can neglect them entierely. 

We do not take into account CP-violating HDOs. These operators are constrained 
by observables such as electric dipole moments. If we choosed to include these CP-odd 
HDOs in our analysis, we would also need to consider the whole set of data sensitive to CP 
violation. Although there is no fundamental problem with such extended analysis, that is 
beyond the scope of the present work. Moreover, the effects induced by CP-violating HDOs 
are often subleading with respect to the effects of CP-even operators, unless the latter are 
sufficiently suppressed. This is the case for Higgs decays, because CP-violating amplitudes 
do not interfere with SM amplitudes, whereas CP-conserving amplitudes do interfere with 
SM amplitudes [19]. 

We could also rigorously take into account the running of the HDO coefficients ai 
from the scale A to the low scale {rrih or ^/s, depending on the process considered), see 
for example [20]. However, the consequences of this running are rather mild so we will 
neglect them in this study. Furthermore, the strong effect of the operator OwB on the /177 
vertex found in [20] requires large enhancement of awB with respect to s^ aww + c^ OiBB — 
\swCwawB (the coefficient of hF^^F^^., see Eq. (2.15)). Given our choice of basis, this 
cannot be explained by a relative loop factor as none of these operators receive contributions 
at tree-level. In the absence of large hierarchies in the couplings of New Physics states, we 



^In fact this is the way how contributions to S and T can arise in theories with new spin-1 states, such 
as in warped extra dimensions [17, 18]. 



conclude that operator mixing does not lead to a large enhancement of the /i — )• 77 rate. 
This result has been shown to hold, independently of the basis of operators [21] (but can 
be somewhat hidden if a basis that does not fulfill our above requirements is chosen). 

2.1 Effective Lagrangian 

In this section we will present the effect of the HDOs on the SM tree-level couplings. 
Loops involving SM particles are considered in Sec. 2.2. We define the physical Higgs field 
h as 

and parametrize the couplings of h to gauge bosons and fermions as ^ 

47 = Xzh {Z^f + Xwh W+W- + Y,Xfh hfR . (2.8) 

/ 

The SM tree-level predictions for these quantities are given in terms of the SM input 
parameters g, v and s^ = g'"^ /{g^ + g'"^)'- 

where the quantities with a tilde are the ones that appear in the SM part of the Lagrangian. 
For instance, g and g' are the couplings appearing in the covariant derivatives. However, 
these couplings do not take the same values as in the SM, since there are corrections from 
HDOs. There are distinct effects, as follows (see Ref. [22] for an analogous discussion on 
fermion couplings). 

• Operators such as 0^2 correct directly the tree-level SM vertices. 

• Some operators {e.g. 0^)2, Oww) modify the kinetic terms of Higgs and gauge fields 
and thus indirectly lead to the rescaling of some couplings. 

• Finally, there can be indirect effects from input parameters. They are taken to be 
the fine-structure constant a, the Z boson mass mz and the Fermi constant Gp, as 
well as the physical fermion masses ruf and the strong coupling constant Us- These 
quantities receive corrections from HDOs but must be held fixed in the analysis. 
Yet, this causes the SM parameters g, v and s^ to become functions of the HDO 
coefficients. 

The last point is sometimes not taken into account in the literature. Let us focus on 
it and define the quantities f , g and Sw via 

Ana = slg\ m| = ^, Gf = ^-,. (2.10) 

^For the hVV couplings with different tensor structure see below. 



These quantities can be viewed as the "famihar" numbers from the SM {e.g. v = 246 GeV). 
Like the input parameters they stay fixed in our analysis. On the other hand, the param- 
eters g, Sw and v are the gauge coupUngs appearing in the covariant derivatives and the 
vacuum expectation value (vev) of the Higgs field, and must be expressed in terms of the 
HDO coefficients. We present the details of this procedure in Appendix A. Taking into 
account all the above effects, we obtain 

Xz = az^, Xw = aw -, Xj = -Cf ^ , 2.11 

V V V 

where mj and rriw are the physical masses. In particular, mw is given by ^ 



m,, 



''^'^l-..^^ + #^V (2.12) 



In the last row we have used Eqs. (5.13) and (5.14) in order to compare our derivation of 
m\Y with the one in [22]. The SM prediction of m\Y is thus only corrected by the oblique 
parameters. In this parametrization, the rescaling factors az, aw and c/ are given by 



az 



aw 



1+ 1 


(I 


1 




1+ 1 


(I 


1 

- -ao 


1 , \ 1-2 

4"^V A^' 


1- 1 


(i"'^^ 


1 

- -ao 


«^ A2- 



c/ = 1- ( -a^2 -|aD -a/ 1 ^. (2.13) 

As a nontrivial consistency check, note that the vector anomalous couplings are rescaled in 
a custodially symmetric way [az = aw) once the custodial-symmetry violating operator 
0^2 is turned off. 

To conclude this subsection we compute the direct tree-level HDO contribution to the 
tensor couplings, 

C'r = C7 ^ {Fi^uf + (:gh (G^.)2 + Cz7 h F^,Z^, + (z h {Z^,f + (w h W+W-, , (2.14) 

which are all zero in the SM at tree-level. One finds 



2 2 1 \ ^ A '^ 

-,-1 — \ ^w <^WW + Cw Ot-BB - IzSwCw awB j T2 ' Cg = "GG J^ ' (2-15) 

CZ7 = ( 'i^CwSyj awW - '^CyjSw aBB - 7^(4 - ^^) OiWB j -^2 ) (2-16) 

1 \ V _ „ V 



Cz — ( c^ oiww + s^ oiBB + -CwSyj oiwB 1 T2 ' CvK — 2 aww t^ • (2-17) 

The first two quantities constitute important corrections to the production and decay of the 
Higgs boson. The last two corrections modify the tensorial structure of the SM Higgs~weak 
bosons coupling in a non-trivial way, which is discussed in detail in Sec. 4. 



^Unlike niz and nif, which are input parameters, the W mass is a prediction in terms of input parameters 
and HDO coefficients. 



2.2 Standard Model loop-induced HDOs 

In this section we compute the Standard Model loop-induced operators relevant for 
Higgs physics. These operators contain indirect modifications due to couplings modified 
by the HDOs considered in the previous subsection. We want to make sure that we do not 
double-count possible New Physics contribution to the Higgs couplings. In order to have a 
well-defined HDO framework at loop-level, we should consider that the HDOs we present in 
Eqs. (2.1)-(2.5) are generated exclusively through New Physics states at leading order, and 
enclose higher-order SM corrections only from irreducible loops. ^ Hence, the modified SM 
loops are not included in the tree-level contributions computed in the previous subsection. 
Our strategy is thus to compute the one-loop corrections to C^^"^^ using the couplings shown 
in Eq. (2.13). 

The one-loop Lagrangian is parametrized as 

^i-ioop ^ ^^ ^ ^p^^^2 ^^^^ (g^^)2 ^ Xz^hF^^Z^, . (2.18) 

Let us decompose these couplings according to the particle in the loop, Aj = ^^^ A^. We 
find 6 

A = aw \ = o T^^ Ay{Tw), (^-19) 



(2.20) 



,W _ ^W,SM _ e aw ,-i ( „ r,2 






2 TW 



Bzy{Tw,Kw)] , (2.21) 



./ _ . ,/,SM_ e2 c.^^e}iTf^-2efsl 



^U = -Cf XT, = I^ V ^/ ,,e. {Bz.irj, .j) - Az,{Tf, .;)) . (2.22) 

where N'i and ej are the number of colors and the fraction of electric charge of the fermion 
running in the loop, respectively. We define Tj = Amf/mf^, Ki = Amf/m^ . The form 
factors Ai,B are given in Appendix B. They are defined so that in the decoupling limit, 
Af^y — )• 1 when r — >■ cx), and Az-y —5- 1, Bz-y — )• when t,k ^ oo. 

2.3 Trilinear gauge boson vertices 

The higher dimensional operators that we are considering also affect charged triple 
gauge boson vertices (TGV) . In the parametrization of Ref . [23] , 

Ctgy = -i e K^F^u W-W+ -igc^ kzZ^, W'W^ -igc^ gf [W+W- - W-,W+] Z^ , 

(2.23) 

^This last point is important for NLO QCD corrections, see Sec. 5. 

®Note that in Eqs. (2.19)-(2.22) only the quantites with a tilde appear. Besides the modified Higgs 
couplings, the HDOs we consider only affect the couplings of the fermions to the W and Z bosons, precisely 
via the oblique parameters 5* and T. The latter would in fact only show up in A^ . However these corrections 
are subleading and rather small (few percents at most), so that it is safe to neglect them. 



the deviations from the Standard Model can be expressed in terms of the HDO coefficients 
as follows: 



1 + 



awB v^ 



2 

A2' 



'" = ' " (2044-4)"'^" + ^(4^^"'"^ + ""0 ^ ' ^'-''^ 

where again some indirect effects from fixing input parameters were taken into account. 
Gauge invariance implies the relation kz = gf — {k^ — l)t'^ and one can check that it is 
indeed fulfilled. We then choose K-y and gf as independent couplings. 

3 Data treatment 

We exploit the results from Higgs searches at the LHC and at Tevatron as well as 
electroweak precision observables and trilinear gauge couplings. The results from Higgs 
searches are given in terms of signal strengths fi{X,Y), the ratio of the observed rate for 
some process X — )■ /i — )• y relative to the prediction for the SM Higgs. An experimental 
channel is defined by its final state (77, ZZ, Z7, WW, bb, tt) and is often divided into 
subchannels having different sensitivity to the various production processes. The accessible 
production mechanisms at the LHC are i) gluon-gluon fusion (ggF), ii) vector boson fusion 
(VBF), in) associated production with an electroweak gauge boson V = W, Z (VH), and 
iv) associated production with a ti pair (ttH) . 

As higher dimensional operators modify not only the Higgs decays but also its pro- 
duction (see Section 5), care has to be taken in extracting the information given by 
the experiments. In particular, when available, we use the results given in the plane 
(/i(ggF + ttH, y),/i(VBF + VH, y)), which (partly) account for the correlations between 
the subchannels. 

The values for the signal strengths in the various (sub)channels as reported by the 
experiments and used in this analysis, together with the estimated decompositions into 
production channels are given in Tables 1-3. Some of the decompositions into production 
channels are taken from [24] . In case of missing information, we take the relative ratios of 
production cross sections for a SM Higgs as a reasonable approximation. To this end, we 
use the latest predictions of the cross sections at the LHC [25] and at Tevatron [26]. 

In both ATLAS and CMS, the Higgs mass is estimated from the two "high- resolution" 
channels: ZZ and 77. In our analysis, the Higgs mass is set to nih = 125.5 GeV (close to the 
combined mass measurement from the two experiments) since it is not yet possible the take 
it as a nuisance parameter without losing the correlations between production channels. 
We consider experimental measurements of the signal strengths as close as possibe to this 
value. 

We take into account the electroweak precision observables using the Peskin-Takeuchi 
S and T parameters [45, 46]. Beyond S and T, the W and Y parameters [47] should be used 



Channel 



Signal strength fi 



mn (GeV) 



/i ^ 77 (4.8 fb"^ at 7 TeV + 20.7 fb 



-1 



Production mode 
ggF VBF WH ZH 

at 8 TeV) [27, 28] 



ttH 



/i(ggF + ttH,77) 
^(VBF + VH, 77) 



1.60 ±0.41 
1.94 ±0.82 



125.5 
125.5 



100% 



60% 26% 14% 



h^ ZZ (4.6 fb-i at 7 TeV + 20.7 fb"^ at 8 TeV) [28, 29] 



/i(ggF±ttH,ZZ) 
/i(VBF ± VH, ZZ) 



1.51 ±0.52 
1.99 ±2.12 



125.5 
125.5 



100% 



60% 26% 14% 



h -^ WW (4.6 fb-i at 7 TeV + 20.7 fb^^ at 8 TeV) [30, 31] 



/i(ggF ± ttH, WW) 
/i(VBF ± VH, WW) 



0.79 ±0.35 
1.71 ±0.76 



125.5 
125.5 



100% 



26% 14% 



h^bb (4.7 fb-i at 7 TeV ± 13.0 fb"! at 8 TeV) [28, 32] 



VHtag 



-0.39 ±1.02 



125.5 



64% 36% 



h-^TT (4.6 fb-i at 7 TeV ± 13.0 ih'^ at 8 TeV) [28] 



//(ggF ± ttH, tt) 
//(VBF ± VH, tt) 



2.31 ±1.61 
-0.20 ±1.06 



125.5 
125.5 



100% 



26% 14% 



Table 1: ATLAS results, as employed in this analysis. The following correlations are included in 
the fit: p-y^ = -0.27, pzz = -0.50, pww = -0.18, prr = -0.49. 



in the HDO framework. However we find that constraints arising from these parameters 
are by far subleading with respect to our other constraints. Experimental values of S and 
T are taken from the latest electroweak fit of the SM done by the Gfitter Group [48]: 
5 = 0.05 ± 0.09 and T = 0.08 ± 0.07 with a correlation coefficient of 0.91. Regarding 
constraints on TGV, we take into account the LEP measurements [49]: 



^ _ n 070+0.044 

nZ — n QO/1+0.022 



(3.1) 



The global likelihood function is defined as the product of the likelihoods associated 
to the various observables, 

L = -^Higgs X Ls^T X -^TGV ! (3-2) 

where -^Higgs is the product of the likelihoods associated to each of the experimental 
(sub)categories, including available correlations. The likelihood associated to the mea- 
surement of an observable O, given as a central value O and a symmetric uncertainty a, is 
modeled by a normal law. 



Lq oc e 



^{0-Oy/2a'2 



(3.3) 



When uncertainties are asymmetric, we use the positive error bar if (0 — 0) > 0, whereas 
we use the negative error bar if (O — O) < 0. Finally, the CMS bound on the decay channel 
h — )• Z7 is implemented as a step function. 



^MZ7 



1 if fiz-y < 9.3 , 
otherwise . 



(3.4) 



Channel 


Signal strength /x 


m/, (GeV) 


Production i 
ggF VBF WH 


node 
ZH 


ttH 


/i ^ 77 (5.1 fb-i at 7 TeV ± 19.6 fh'^ at 8 TeV) [33] 


/i(ggF ± ttH, 77) 
/i(VBF ± VH, 77) 


0.49 ±0.39 
1.65 ±0.87 


125 
125 


100% 

60% 26% 


14% 


— 


h^ZZ (5.1 fb"^ at 7 TeV ± 19.6 fb"^ at 8 TeV) [34] 


^(ggF ± ttH, ZZ) 
/i(VBF ± VH, ZZ) 


0.99 ±0.46 
1.05 ±2.38 


125.8 
125.8 


100% 

60% 26% 


14% 


— 


h^W\ 


M (up to 4.9 fb-i at 7 TeV ± 19.5 Ib'^ at 8 TeV) [35-38] 


0/1 jet 

VBF tag 

VHtag 


0.76 ±0.21 

-0 31+2-22 


125 

125.8 
125.8 


97% 3% 
17% 83% 

64% 


36% 


- 


h^hh (up to 5.0 fb-i at 7 TeV ± 12.1 fb"! at 8 TeV) [38-40] 


z(r^+)h 

Z{yv)\i 
W(£iy)h 
ttH tag 


1 ^^+1-20 
1 70+1-11 

0.691°!^ 
-0.80+?i° 


125 
125 
125 

125.8 


100% 


100% 
100% 


100% 


h^TT (4.9 fb-i at 7 TeV ± 19.4 ft-^ at 8 TeV) [41] 


0/1 jet 

VBF tag 

VHtag 


7fi+0-49 

1 40+0-60 
-L-40_o.57 

0.771^:^8 


125 
125 
125 


76% 16% 4% 
19% 81% 

64% 


3% 
36% 


1% 


/i ^ Z7 (5.0 fb-i at 7 TeV ± 19.6 fb^^ at 8 TeV) [42] 


Inclusive 


< 9.3 at 95% CL 


125.5 


87% 7% 3% 


2% 


1% 



Table 2: CMS results, as employed in this analysis. The following correlations are included in the 
fit: p^-y — —0.50, pzz = —0.73. 



Channel 


Signal strength fi 


ruh (GeV) 


ggF 


Production 
VBF WH 


mode 
ZH 


ttH 


/i ^ 77 [43; 




Combined 


6 izt+3-25 
0--L4-3.19 


125 


78% 


5% 


11% 


6% 


- 


h -^ WW [43] 


Combined 


omtlil 


125 


78% 


5% 


11% 


6% 


- 


h^bb [44] 


VHtag 


1 56+0-'^2 
i.OD_073 


125 


- 


- 


62% 


38% 


- 



Table 3: Tevatron results for up to 10 fb at \/s — 1.96 TeV, as employed in this analysis. 

We will now derive the deviations induced by the HDOs to the observables presented 
in Sec. 3. We first discuss the particular treatment of tensorial couplings. All formulas are 
given in the following section. 
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4 On weak bosons tensorial couplings 

Because of electroweak symmetry breaking, the W, Z = V bosons generally couple to 
the Higgs through two different Lorentz structures. The coupling can be vectorial, oc g^'^, 
or it can be tensorial with a vertex oc {gf^'^ — ^i^^^^ where qi, q2 are the momenta of the two 
gauge bosons. The leading SM couplings \w, ^z given in Eqs. (2.8) and (2.9) are vectorial. 
Tensorial couplings are generated only at one- loop and are 0{a) ~ 10~^. 

Once HDOs are taken into account, the relative importance of the vectorial and tenso- 
rial terms is modified. On one hand vectorial couplings are rescaled by the coefficients a\y^z- 
On the other hand new tensorial contributions (w^ Cz are generated following Eq. (2.17). 
The amplitude associated to a hVV vertex (with the V^s possibly off-shell) is in general 



M{hVV) 



Ai,A2 _ gA«(*)g'^(*) 
Ai A2 



^av\f^g^'' 



i2Cvqi-q2 



flU 



91-92 



(4.1) 



where M'^'^ and A^^'^ are the longitudinal and transverse helicities amplitudes, respec- 
tively. Interferences among helicity amplitudes then determine angular distributions (see 
e.g. [50]). In this work, we consider that the SM contribution to the tensorial coupling 
is small with respect to the one induced by New Physics. The relative magnitude of the 
longitudinal and transverse amplitudes in case of a vectorial coupling is given by 



A^r 



Mt^^ 



mt 



qI 



92 



2kilk2| 



(4.2) 



while it is the inverse in case of a tensorial coupling, 

2kilk2| 



n 



Mr 



M 



±,± 



\mt 



^l 



'/2I 



(4.3) 



The two vector bosons can be off-shell in the above expression, while the Higgs is on-shell. 

As r^ y^ Tt, the two Lorentz structures imply generally different angular distributions. 
Moreover, even for unpolarized processes, the energy dependence in Eq. (4.1) is different 
for both contributions, such that also energy distributions are modified. Because of this 
different energy dependence, kinematic cuts prepared for the SM are generally unadapted to 
such a non-trivial modification. That is, in Eq. (5.1), esM+HDO ¥" ^SM- The consequences 
may be an incorrect estimation of the signal strength and of the Higgs mass. To perform 
an exact analysis, one should redo the fits to LHC data taking into account the modified 
Lorentz structure in the expected signal. Such work is clearly beyond the scope of our 
present study. Instead we will show that under reasonable approximations we can use 
esM+HDO = ^SM in the present analysis. 

There are three processes sensitive to the Cv tensorial couplings in the context of 
the searches for the Higgs boson at around 125 GeV: the leading decay to weak bosons 
h —7- VV* , and the VBF and VH production modes. We now discuss how we treat these 
three tensorial contributions. 
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4.1 h -^ VV* 

In the case of a light Higgs boson, the leading decay occurs with one of the V off 
the mass shell. The weak bosons then decay into fermions. For massless fermions, the 
kinematic bounds on the on-shell boson energy Ey are my < Ey < {fn\ + 'rny)/'2,mh in 
the rest frame of the Higgs. Because of the V* propagator, the lower bound Ey = rny is 
favorized, implying that both weak bosons are preferentially produced at rest. Longitudinal 
and transverse amplitudes are then equally populated, r„ = 1. Therefore, one has r^ = 1 
as well, such that one can see qualitatively that a tensorial contribution cannot radically 
modify angular distributions. This is confirmed with the exact angular and invariant 
mass distributions among leptons induced by pure vectorial and pure tensorial couplings 
[51, 52]. ^ In our study, the tensorial contributions are constrained to be subleading with 
respect to the vectorial contributions, such that the deviations induced on angular and 
invariant mass distributions can be easily be smaller than the current statistical uncertainty. 
In addition, they could also be misidentified with the background. For example, in /i — )• 
VV* , the distribution of the most discriminant observable, "lepton-opposite Z momentum 
angle" , is very similar to the distribution of the irreducible background qq — )• ZZ* (see 
Fig. 3 in [53]). 

Following what discussed above, we can reasonably assume that angular and invariant 
mass distributions are not affected by the presence of tensorial couplings given the current 
level of precision. Polarization of the on-shell V can thus be averaged, and we are left with 
a matrix element scaling as 

\M\^ = \My+Mt\^ cc \ayXf^ - 2Cvqi.q2f , (4.4) 

where qi, q2 are the momenta of the two vector bosons. In the Higgs rest frame, one has 
qi.q2 = rrihEy — rriy, which is bounded as 

2 2 

/ N 'Tift '^V / A r\ 

my{mh-my) <qi.q2< -. (4.5) 

The exact tensorial contributions to the total decay widths are given in App. C. We intro- 
duce the dimensionless positive quantity 

vvv = qi-q2/ml, (4.6) 

with V = W,Z. Defining 

JuyyM.MtdPS 2 ^ f iy^yy\Mt\^dPS 

^"^^^= jM.M^dPS ' ^"^^^= J\Mt\^dPS ' ^^-^^ 

the vector-tensor interference term will be oc Cvi^vv) and the pure tensor contribution 
will be oc \Cv\'^{i^vy)- For rrih = 125.5 GeV, mz = 91 GeV, mw = 80GeV, one gets 
{uzz) = 0.2209, (4z)^ = 0.2211, {vww) = 0.2653, {i^ww)^ = 0-2659. In the following we 
will make the approximation {I'yy) ~ (uyy)'^. 



^ Overall, the situation is much less striking than for a CP-violating contribution, which forbids the 
decay to the longitudinal polarization state. 
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4.2 VBF production mode 

For the VBF process, both ATLAS and CMS apply hard cuts on the outgoing jets 
rapidities and their difference. The rapidity distributions of the two jets are similar in 
presence of a tensorial coupling, just like in the decay into two photons or in the production 
via gluon-gluon fusion, such that one can assume that cut efficiency is the same. The crucial 
change lies in the azimuthal angle (pjj between the two tagging jets (see e.g. [50] and 
references therein). Indeed, both weak bosons are space-like, with virtualities considerably 
smaller than m|. Such values are favorized to balance the space-like V and the outgoing 
jets virtualities. As a result one has typically r„ 3> 1, rt <C 1 i.e. vectorial and tensorial 
amplitudes are mostly longitudinal and transverse, respectively. Consequently, the cpjj 
distribution is almost flat for a pure vectorial coupling, and strongly peaked at tt/2 for a 
pure tensorial coupling. For a large enough HDO contribution to the tensorial coupling, 
an anomalous (j)jj distribution could thus be observed. However, this variable is not used 
for the selection of the events in the experimental analyses we consider. Therefore, the 
selection efficiencies are also suitable in the case of large tensorial contributions, and one 
has £sM = esM+HDO- One can average over the polarizations, and the squared amplitude 
is then simply rescaled by a factor \av^]/BF ~ '^CvQi-(l2\ ■ 

We still have to determine the magnitude of the tensorial contribution. In this process, 
the scalar product of the weak boson momenta qi.q2 is related to the incoming and outgoing 
quarks as qi.q2 = "t-I/^ + Pi-P3 + P2-Pa- The outgoing quarks are highly energetic with 
respect to the amount oi pT they receive from the V fusion, such that one has |pi| ~ IpsI 
and I p2 1 — I P4 1 . In terms of the pT and rapidities of the outgoing quarks we have then 

qi-q2 = ^ + \PT,zr — 2 — + \pta\ — 2 — • ^^-^^ 

Without the tensorial contribution, the pT distribution peaks typically at values smaller 
than mv- The tails of the pT distributions drop quickly for higher energies [55], with 
typically one jet at a time getting a large pT [56]. One can thus assume qi.q2 ~ ""T-h/^ ^° ^ 
good approximation. Once the tensorial coupling is taken into account, a deviation from 
the expected SM distributions might be present in the high-p^ tails, as qi.q2 is enhanced 
at large pT- However, as long as one counts the total number of events, i.e. the integral of 
the distribution, this enhancement of qi.q2 has a small weight and can be safely neglected. 
Finally, defining the dimensionless positive quantity 

vvBF = qi-q2/'m\, (4.9) 

with V = W, Z, we have thus ^vbf ~ 1/2 after phase space integration. 

4.3 VH production mode 

In the case of the associated production with an electroweak gauge boson, the scalar 
product of the momenta of the weak bosons is given by 



s + m?/ — m? 
qi-q2 = ^ ^, (4.10) 
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where ^/s is the partonic center-of-mass energy, which is typically s/niy = 0(100) at the 
LHC. Therefore, contrary to the two other processes, the product qi.q2 is large because it 
contains the partonic center of mass energy ^/s. The tensorial contribution is then sub- 
stantially enhanced in this process. Besides, we have r^ 7^ r^, as r^ = r^ ~ y/sjlmy^ such 
that the angular distributions may in principle be substantially modified by the presence 
of the tensorial coupling. 

However, it turns out that for both polar and angular distributions, the angular effects 
can be neglected. We refer to [57] and references therein for the expressions. Although 
results are given for e^e~ collisions, they can be trivially generalized in the case of the LHC. 
For the distribution of the polar angle of the vector boson in the laboratory frame, it is 
the longitudinal component of V which enters mainly, such that the tensorial contribution 
to the distribution is suppressed by an additional factor 0{my/s). For the azimuthal 
distributions, the tensorial contributions can be sizeable, but the whole distribution tends 
to be flat for s >> my, with non-flat terms suppressed by powers of mv/\/^. As a result, 
although various pieces of angular information are used in event selection for this mode of 
production, we can safely neglect the angular effects of the tensorial coupling. 

Concerning the magnitude of the tensorial contribution, it appears that it reduces to 
a simple rescaling oc Xy^ + \2C,v'my in the limit s » my. The rescaling is exact up to 
a subleading term 0{12'my/s) « 0.1. To include the subleading s-dependent terms, an 
integration over the partonic density functions would be necessary. 

5 Deviations caused by New Physics 

5.1 Higgs signal strengths 

Theoretical signal strengths for Higgs searches can be expressed as 

uiX Y) - [-(-^^^)^(^^^)^^1sM^HDO .5 ,. 

^^ ' ^" [a{X^h)B{h^Y)e^yU ^^^ 

where B is the branching ratio of the decay and the coefficient e^^ G [0, 1] characterizes the 
efficiency of event selection for a given subcategory. In all generality, efficiencies in the SM 
with and without HDOs are not necessarily the same, i.e. ssm+hdo / ^SMj because kine- 
matic distributions can be modified in a non-trivial way by HDOs. The selection criteria 
calibrated on the SM expectations are then unadapted in such situation and complicates 
the interpretation of the signal strengths. 

However, we have seen in Sec. 4 that one can safely ignore these possibilities of HDOs af- 
fecting the kinematic distributions, given the current precision of the experimental searches. 
It is therefore a good approximation to set esM+HDO = £SM- Thus, for each signal strength, 
one can simply incorporate the contributions coming from the tensorial couplings in the 
rescaling of the Standard Model signal strength. 

The gluon-gluon fusion process is modified both by the tree-level HDO contribution Q 
and the anomalous Higgs-fermion couplings c/. Keeping only the third generation, we get 



(TggF — CTggp 



\t,SM , xfe.SM , > 

CtXg + CbXg + Cg 



-,t,SM I xfe,SM 
Ag + Ag 
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(5.2) 



Vector boson fusion is modified by the anomalous vectorial couplings avF,z and by Cw,z- 
Denoting by Ay^p the effective SM couplings, one has 



C^VBF = CTvBF 



iSM 



aw^W + «^^z - 2 lA/BF mi iCw + Cz) 



\SM I \SM 



(5.3) 



The parameter ^vbf is defined in Sec. 4.2. We take i/vbf = 1/2- The associated production 
with an electroweak gauge boson is modified as 



(TVH = CTvH 



av\f^ + l2Cvm\ 



ASM 



where V = W, Z. Finally, the associated production with a tt pair is rescaled as 

0"ttH 

The decays of the Higgs boson into fermions are modified as 



U |2^SM 



// 



\cf\'r\ 



2-pSM 
// 



The tree-level decays to vector bosons are modified as 



Tvv 



avXf'-2Cvml{uvv) 



\SM 
Ay 



-nSM 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



where the parameter (i^yy), defined in Eq. (4.7), encodes the modification of phase space 
integrals. Loop-induced decays are sensitive to more deviations, 



■77 



Z7 




\W,SM , \t,SM , ,6,SM , .T,SM 
^Zj + ^^7 + ^Z-y + ^Z-y 



(5.8) 

(5.9) 



In such cases, the tensorial couplings can compete with the SM effective couplings. 



5.2 QCD radiative corrections 

Many of the above described processes receive leading radiative corrections from QCD 
loops. For all the tree-level processes, the structure of loop diagrams is not modified by 
the insertion of HDOs, including the tensorial couplings, such that radiative corrections 
factorize up to higher order corrections. It is thus straightforward to take them into 
account, simply using the NLO predictions of a and F . 

The situation is more involved in the case of the loop-induced processes {h — )■ 77, 
h —7- Z7, and gg — )■ h) because this time the tensorial coupling is competing with the 
SM loops. Hence the effects of the C's may be very large in these processes, such that it 
is important to properly take into account the radiative corrections. As stated in Sec. 2, 
the HDOs implicitely contain higher-order corrections from irreducible SM loops. These 
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contributions therefore have to be taken into account for the SM effective couplings and 
not for the C, couphngs. ^ 

The processes /i — t- 77 and h — t- Z7 only receive virtual NLO QCD corrections. For 
h —7- 77, we take into account the exact values of the correction factor to the quark effective 
couplings 



Af^^ = Af^|Lo 1 + -C^(r,)) , (5.10) 



vr 

where the Ch function can be found in [56]. For h — )• Z7, one can take the correction in 
the heavy top limit as a good approximation [56], 



a, 

vr 



Xf^ = X^^^\^oil-^)- (5.11) 



The situation is more subtle for the ggF process, because of the presence of important 
NLO real corrections. Introducing the tensorial coupling leads generally to non-trivial 
modifications of the integrals over parton densities for real emissions. However, in the 
heavy-top limit and neglecting the small bottom quark contribution, the QCD corrections 
to the SM loop and to the tensorial coupling (^g become similar and factorize. Adopting 
this fairly good approximation, the SM effective coupling are rescaled as 

Ar = Af-|.o(l + ^f). (5.12) 

5.3 S and T parameters 

The electroweak precision observables are affected in the presence of the HDOs. At 
tree-level the S and T parameters are related to the HDO coefficients as follows: 

aS = (2 SyjCyj awB + s^ "£> + ^^ "d ) T2 ' (5.13) 

aT = f--a^2 + -a'z)j ^. (5.14) 

Moreover, the SM loops are modified by the HDOs. The T parameter receives new divergent 
contributions from the modified SM couplings az and ay/ in Eq. (2.13). A quadratic 
divergence, 

A^ a'r^2v'^ 
aAT = --— ^^-L 5.15 

arises from custodial breaking [58] . Dropping other terms that are proportional to a^ and 
a'^2 (that already appear at tree-level) the two couplings coincide and we can take the 
result from Ref. [5], 

3e^ / 1 \ v"^ fmh\ 

aAT = - „^ ^ ^ a£,2 - -ao t^ log — - . (5.16) 



32 7r2 4 V 2^;a2^^va; 

Similarly, the S parameter receives corrections due to the modified Higgs coupling az [5], 
hence it is expected to get new contributions proportional to a£)2 and a'jj2- Finally, the 



^We are grateful to M. Spira for enlightening discussion on this subject. 
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tensor couplings C,y can also generate new SM loop contributions which have been given 
in Ref. [59], 

6/ 1 \ V / Tfih \ 6 V / Tflh \ 

"^"^ " 24^ V"^' ^ 2"'^7 A2 ^°^ VX; ^ 2^^"-^^ + aww) j^^og (^— J . (5.17) 

Finally we neglect the contraints coming from the W and Y parameters [47] as they are 
expected to have a small impact on our results. 

6 Bayesian setup and low-A scenario 

6.1 Bayesian inference 

We are working in the framework of Bayesian statistics (see [13] for an introduction). 
In this approach, a probability is interpreted as a measure of the degree of belief about 
a proposition. Our study lies in the domain Bayesian inference, which is based on the 
relation 

p{e\d,M)o^p{d\e,M)p[e\M), (e.i) 

where 9 = {0i...n} are the parameters of the model A^, and d denotes the experimental data. 
The distribution p{6\d,M.) is the so-called posterior probability density function (PDF), 
p{d\9,M.) = L{0) is the likelihood function enclosing experimental data, and p{9\M.) is the 
prior PDF, which represents our a priori degree of belief on the parameters. The model 
M is in our case the Standard Model extended with higher dimensional operators. The 
likelihood is defined in Sec. 3 (see Eq. (3.2)) and the theoretical expressions for the HDO 
modified signal strengths are given in Sec. 5. The prior PDF is discussed in the next 
subsection. 

The posterior PDF is the core of our results. Integrating the posterior over a subset A 
of the parameter set ^ = {'0, A}, 

p(V'|d,A^)oc f d\p{i;,\\M)L{tP,X), (6.2) 

leads to inference on the parameters ■(/'• 

Also the notion of naturalness and fine-tuning are built-in [14] in the Bayesian ap- 
proach. This is relevant for our study, in which precise ("fine-tuned") cancellations be- 
tween various HDO contributions can happen. Intrinsically, the regions of parameter space 
in which precise cancellations occur have a weak statistical weight, such that they are 
flushed away after integration. The results we will present can thus be considered as 
generic, i.e. free of improbable cancellations. 

We will consider uniform (flat) priors for the quantities 

and demand |/3j| < 1. Moreover, we will flx the cutoff scale to be A = Attv. In the following 
we will justify these choices and argue that it ensures in particular convergence of the HDO 
expansion as well as perturbativity of the UV theory, and minimizes the dependence on 
the choice of the HDO basis. 
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6.2 Priors and low-A scenarios 

The prior distributions associated to our parameters is a key feature of Bayesian infer- 
ence. We follow the "principle of indifference" [60, 61] that maximizes the objectiveness of 
the priors. Once a transformation law 7 = f{6) irrelevant for a given problem is identified, 
this principle let us find the most objective prior by identifying pQ = pr in the relation 
Pe{0)d9 = pr{-f)dj. 

The cutoff scale A is given a logarithmically uniform PDF, 

p(A) « J . (6.4) 

By doing so, all order of magnitudes are given the same probability density. Regarding the 
dimensionless coefficients a, note that the choice of the HDO basis should be irrelevant for 
the conclusions of our study. Given that coefficients in different basis are related through 
linear transformations, the most objective prior to associate to each Oi is the uniform 
PDF, 9 

p{ai) oc 1 . (6.5) 

This choice of prior is well justified, however, one should keep in mind that other possibil- 
ities still exist. 

Let us emphasize that in our general framework, the following hypotheses need to be 
scrutinized. 

• Perturbativity of the HDO expansion, |aj|/A^ < 0{l/v'^), 

• Perturbativity of the couplings expansions in the UV theory, |aj| < 0(167r^), 

• HDO generation by loops, 

• Custodial symmetry. 

In the present work, we investigate scenarios of low-scale New Physics, with values 
of A going up to 0{A7rv). We take custodial symmetry to be an exact symmetry of the 
theory. This forbids the presence of the operators 0^2 and O^. As a consequence, one 
has aw = az = ay and some contributions to the EW precision observables are suppressed 
including the potentially large quadratic divergence in T. Recall that Oww-, (^WB, and 
Obb are all independently custodially symmetric. This generally implies that processes 
involving the W and Z are not identically rescaled, for instance 

CWR , q"ZH .„ „x 

_SM 7^ ^SM • v"-"'' 

"WH "ZH 

Our approach goes therefore beyond the fits involving pure rescalings induced by anomalous 
couplings. 



® Here the indifference principle sets the shape of the PDFs but does not set the bounds. One can see 
that ranges on a's are not conserved from one basis to another. In the scenario of democratic HDOs, this 
issue will be automatically solved, as one relies only on perturbativity of the HDO expansion to set the 
bounds on a's. In the scenario of loop suppressed Off's, one takes advantage of a particular choice of 
basis, so the same argument does not apply in that case. 
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I) Democratic HDOs 


II) Loop-suppressed Off's 


A 


47rz; 


Airv 


Pff 


[-1,1] 


[-1/I67r2,l/I67r2] 


Other /3 


[-1,1] 


[-1,1] 



Table 4: Summary of the setup of the scan in the two scenarios we consider. The f3pp = app v^ / A? 
coefficients (where FF = WW, WB, BB, GO) correspond to the field-strength-Higgs operators, 
fn both cases we take custodial symmetry to be an unbroken symmetry. 

Over this range of A, perturbativity of the HDO expansion is the dominant constraint 
as it requires |aj| < b^ jv^ which automatically implies |aj| < IGvr^ and hence perturbativ- 
ity of the couplings expansions in the UV theory. 

When the HDOs are generated within a perturbative UV theory, none of the field 
strength-Higgs operators Off = Oww,wb,bb,gg (see Eqs. (2.2) and (2.3)) can be gener- 
ated at tree-level. Because of our appropriate choice of basis, these loop-generated HDOs 
are exactly the ones associated with the tensorial couplings Cg,7,Z7- We will therefore distin- 
guish between two scenarios, depending on whether or not the Off^s are loop suppressed 
with respect to the other HDOs. Given that tensorial couplings can play an important role, 
this distinction is particularly crucial. The two scenarios, denoted by I and II, are respec- 
tively dubbed "democratic HDOs" and "loop-suppressed Off^s". The main features are 
summarized in Table 4. We emphasize that the interpretation of A as a true New Physics 
scale also depends at which order the whole set of HDOs is generated. For instance, in the 
i?-parity conserving MSSM, the whole set of HDOs is generated only at one-loop order, 
such that the actual NP scale should be 0(47rA). 

A parameterization particularly adapted to low- A scenarios is as follows. Defining the 
parameters 

Pi = Oiij^, (6.7) 

it follows that the /3's and A are independent, i.e. p(aj. A) = p(A)p(/3j). The /3's prior is 
the uniform PDF over [—1; 1], noted U{/3i). The prior of A is p(A) oc A^""^, where n is the 
number of /3's. In our case, n = 9 is large enough such that this prior is essentially peaked 
at Amax, p(A) ~ 5{A — Ajnax)- We have therefore 

piau A) = 5(A - A„,ax)C/(/3i) . . . C/(/3n) • (6.8) 

This factorization allows us to marginalize over A, and to present our results in terms 
of /3's, which contain all the relevant information. A mild dependence on A will remain 
through loop-level ©(log A) terms in the S and T parameters, that will be discussed below. 
The fact that /3's prior is uniform and spans a constant range is essential to facilitate 
interpretation of the posterior PDFs. The fact that A ~ Amax is also useful, as it renders 
straightforward the evaluation of the few A-dependent terms. 

This parameterization turns out to be convenient in order to extract information about 
HDOs in a scale independent way, up to a mild ©(log Amax) dependence. For example, for a 
given A, one can directly read the values of a's on the /3's plot. Similarly, for given a's, one 
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can deduce the allowed A values from the plots. This parameterization is appropriate at 
low A, up to A = 0{4:TTv). Beyond this scale, the bound from HDO perturbative expansion 
competes with the bound from the perturbative expansion of the couplings. Once the latter 
dominates, the features of factorization no longer hold. 

6.3 The MCMC setup 

We evaluate posterior PDFs by means of a Markov Chain Monte Carlo (MCMC) 
method. The basic idea of a MCMC is setting a random walk in the parameter space such 
that the density of points asymptotically reproduces the posterior PDF. Any marginal- 
isation is then reduced to a summation over the points of the Markov chain. We refer 
to [13, 62] for details on MCMCs and Bayesian inference. Our MCMC method uses the 
Metropolis-Hastings algorithm with a symmetric, Gaussian proposal function. We run re- 
spectively 50 and 15 chains with C(IO^) iterations each for the democratic HDOs case and 
the loop-suppressed OpF^s case. Finally, we check the convergence of our chains using an 
improved Gelman and Rubin test with multiple chains [63]. The first 10 iterations are 
discarded (burn-in). 

7 Inference on HDOs 

In this section we present and analyze the posterior PDFs arising in our scenarios 
of democratic HDOs and loop-suppressed OpF^s, denoted by I and II, respectively. Our 
results will be shown in terms of the /3j = ajW^/A^ parameters, which encode information 
about the fundamental parameters. Recall that the /3's prior PDF is uniform, and that 
the /3 PDFs we show are valid for any value of the cutoff scale A < 0{4:Ttv), up to a mild 
log A dependence. Moreover, this parameterization sets A ~ Amax- The posterior PDF we 
present is computed for A = Airv ~ 3TeV. For smaller A, we expect the oc log A constraints 
from AS and AT to mildly relax. We will comment below on this effect. 

We will also discuss deviations from the SM cross sections and decay widths, defining 

"X ^ h^Y ^ h 

where X = ggF, VBF, WH, ZH,ttH, and Y = j-f, ZZ, Zj, WW, bb, tt. Note that 
the observables are the signal strengths ft{X,Y) rather than the individual Rx and Ry- 
Furthermore, the total width of the Higgs boson is about 4 MeV in the SM and cannot be 
probed directly currently at the LHC. The signal strengths, associated with a production 
mechanism X and a decay Y, can be expressed as 

f,iX,Y) = Rx-^. (7.2) 

We present one-dimensional PDFs of the fundamental parameters /3j for both scenarios 
in Fig. 1. Moreover, in Table 5 we report the 68% and 95% Bayesian credible intervals 
(BCIs) for these quantities. We also present the BCIs for the other, dependent quantities, 
i.e. the anomalous couplings oy and Cf, the tensorial couplings d, and the various i?'s. 
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Figure 1: Posterior PDFs of the 9 fundamental parameters, (3i = a^w^/A^, in scenario I (black) 
and scenario II (red). 



One can first remark that all of our HDO coefficients except (3t and Ph are constrained 
enough to stay within the bound |/3j| < 1, as required for the convergence of the HDO 
expansion. Furthermore, the Pff = Pww,wb,bb,gg coefficients are 0(0.01) in both sce- 
narios. Pd and PwB are strongly correlated in both scenarios as they appear in the S 
parameter at tree-level, see Eq. (5.13) (we recall that we fix a'jy = a'^2 = in order to pre- 
serve custodial symmetry). We thus have 2 Cw f^wB ~ —Sw Pd as can be seen in Fig. 2. The 
TGV observables also involve (3d and /3vkb (see Eq. (2.24)), and thus provide an indepen- 
dent constraint on f^jj (or equivalently Pwb)- The slight deficit in K-y and gf as measured 
by LEP, see Eq. (3.1), tend to favors positive (negative) /3wb (Pd)- Finally, note that in 
scenario II the PDF of PwB is limited to the [— l/lGvr'^, l/lGvr^] range since we consider 
that the operator OwB is loop-suppressed. This in turn fixes the allowed range for Pd. 

The (3^2 coefficient is allowed to deviate significantly from as it only appears in 
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Figure 2: Posterior PDFs of /3iys versus /3d in scenario I (left) and scenario II (right). The red 
and blue regions correspond to the 68% and 95% Bayesian credible regions (BCRs). The green star 
indicates the maximum of our posterior PDF. 

loop contributions to S and T and in ay. The probability of having /3£)2 > is 94% 
(90%) in scenario I (II) and comes from T, as well as VBF and VH production modes and 
h —7- VV decays. A value for ay > 1 leads to a positive contribution to T, as well as an 
enhancement of the VBF and VH production processes, the h — t- VV* decays, and also to 
the loop-induced decay rates, /i — >• 77 and h — )• Z7. 

Pww and Pbb are mainly constrained by the searches for /i — )• 77 and h — )■ Zj as 
they contribute to the tensorial couplings C7 and (z-y, see Eq. (2.15) and (2.16). Given the 
large allowed range for /3vkB) a cancellation has to occur with Pww and fSsB in order to 
achieve a /i — )• 77 rate compatible with experiment. As a result, negative values of Pww 
and /3bb are favored. Moreover, Pww is also constrained from VH production processes via 
the quantities (y In contrast (3bb and /3\yb play no role for these measurements, as they 
are more strongly constrained by the other effects mentioned above. The contributions of 
(3ww and Pbb to S are up to 0(0.03) and do not impact the PDFs. This effect is even 
smaller in scenario II, and is also smaller if we take A < Attv due to the log(m/i/A) factor 
in Eq. (5.17). We note that in scenario II the PDFs for /3vKVF and /3wb can easily reach 
the bounds set by the priors, while Pbb is more strongly constrained by the data. This is 
due to the fact that Pbb enters in (^^ with a coefficient roughly four times larger than the 
other two. 

Finally, the Yukawa corrections parametrized by /?/ (/ = t, b, r) are much less con- 
strained as they only contribute to the rescaling factors c/, but account for most of the 
deviations of c/ from 1, such that we often have \(3f\ » |/3d/4| and thus Cf ^ 1 + (3f. It 
is worth noting that f3t has a fairly large probability of being close to —1, which leads to 
small or vanishing q. The posterior PDF of q is shown on the left pannel of Fig. 3. 

In such case, one may wonder whether or not the preference for small q is due to a 
volume effect caused by the process of marginalization. To this end, we display on the 
right panel of Fig. 3 the profile likelihood for the parameter q, i.e. the likelihood for given 
cj, maximized over all the other parameters. We conclude that in both scenarios, the 
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scenario I 


scenario II 




68% BCI 


95% BCI 


68% BCI 


95% BCI 


Pd X 103 


[10, 120] 


[-50, 180] 


[-6,23] 


[-19,26] 


/3o2 X 10^ 


[70, 350] 


[-50, 480] 


[40, 290] 


[-90, 400] 


A X 10^ 


[-1000,110] 


[-1000,610] 


[-930, 10] 


[-1000,590] 


15b X 103 


[-10,530] 


[-220, 930] 


[-110,500] 


[-280,860] 


Pr X 10^ 


[-170,300] 


[-420,510] 


[-190,270] 


[-450,510] 


/3gg X 103 


[-3.2,8.0] 


[-4.0,9.6] 


[-3.3,0.6] 


[-4.2,2.7] 


/3w-w- X 103 


[-19,7] 


[-30, 18] 


[-5.6,2.3] 


[-6.0,5.6] 


^W'i? X 103 


[-32,1] 


[-49, 13] 


[-6.0,1.6] 


[-6.3,5.3] 


/3iJij X 103 


[-12,0] 


[-17,4] 


[-1.7,1.6] 


[-2.9,3.0] 


av 


[1.02,1.15] 


[0.96,1.21] 


[1.02,1.14] 


[0.96,1.20] 


Ct 


[0.05,1.14] 


[0.03,1.63] 


[0.06,1.01] 


[0.04,1.60] 


Cb 


[0.90,1.54] 


[0.79,1.96] 


[0.89,1.50] 


[0.72,1.86] 


Ct 


[0.84,1.31] 


[0.58,1.53] 


[0.81,1.27] 


[0.55,1.51] 


CgVX 103 


[-3.2,8.0] 


[-4.0,9.6] 


[-3.3,0.6] 


[-4.2,2.7] 


C^VX 103 


[-5.5,0.5] 


[-6.1,0.9] 


[-0.33,0.46] 


[-0.69,0.86] 


CzyVX 103 


[-13,18] 


[-18,30] 


[-4.9,4.4] 


[-7.6,7.9] 


Cz f^ X 103 


[-20,2] 


[-31,11] 


[-3.4,2.3] 


[-5.1,4.4] 


CvK -(^ X 103 


[-39, 15] 


[-59,37] 


[-11,5] 


[-12,11] 


i?ggF 


[0.6,1.3] 


[0.5,2.0] 


[0.6,1.3] 


[0.4,2.0] 


-RVBF 


[1.0,1.4] 


[0.9,1.6] 


[1.0,1.3] 


[0.9,1.4] 


i?WH 


[0.7,1.3] 


[0.5,1.7] 


[1.0,1.3] 


[0.9,1.4] 


i?ZH 


[0.7,1.2] 


[0.5,1.5] 


[1.0,1.3] 


[0.9,1.4] 


^ttH 


[0.02, 1.0] 


[0.02,2.6] 


[0,0.9] 


[0,2.5] 


^77 


[1.1,1.9] 


[0.8,2.5] 


[1.1,1.8] 


[0.8,2.3] 


Rzj 


[0,5.2] 


[0,12.0] 


[0,2.2] 


[0,4.3] 


Rzz 


[1.0,1.3] 


[0.9,1.5] 


[1.0,1.3] 


[0.9,1.4] 


Rww 


[1.0,1.3] 


[0.9,1.5] 


[1.0,1.3] 


[0.9,1.4] 


Rbb 


[0.7,2.2] 


[0.5,3.6] 


[0.7,2.1] 


[0.4,3.3] 


Rtt 


[0.6,1.6] 


[0.3,2.2] 


[0.6,1.5] 


[0.2,2.1] 


-Rwidth 


[0.8,1.9] 


[0.7,2.7] 


[0.8,1.8] 


[0.6,2.5] 



Table 5: 68% and 95% Bayesian credible intervals (BCIs) for the democratic HDOs case (scenario 
I) and for the loop-suppressed Opp^s case (scenario II). 

preference for small q originates from the likelihood and not from a volume effect. ^^ 

The shapes of the PDF and profile likelihood for ct in Fig. 3 are in fact a direct 
consequence of the signal strength measurement /i(ttH, bb) by CMS [40], see Table 2. Notice 
that the latter is so far the only analysis sensitive to the ttH production mode. In spite 
of its large error, the low central value drives q efficiently to small values because of the 



'^We have also checked that no volume effects appear for any of the other posterior PDFs either. 
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Figure 3: On the left, posterior PDF of ct in scenario I (black) and scenario II (red). On the right, 
profile likelihood along the ct axis in scenario I and scenario II (same color code). 
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Figure 4: Posterior PDF of /3gg versus /3t in scenario I (left) and scenario II (right). Color code 
as in Fig. 2. 

relation i?ttH = c^ . Although small q decreases (increases) the value of i^ggF (-R77)) these 
changes can be compensated for without decreasing the likelihood. In the case q w 0, the 
gluon-gluon fusion (ggF) process is mainly driven by the tensorial coupling Cg = I^gg/v. 
We show in Fig. 4 the correlation between ficG ^-nd /3i, which is needed to reproduce the 
observed ggF rate. For the decay h — )• 77, we observe an increased rate R-y^f > 1, which 
can be seen in Fig. 5. Indeed, in the SM the /i — t- 77 process is dominated by the W loops, 
and there is a destructive interference between the t and W contributions. Therefore, the 
suppression of cj helps increasing R^^. To better understand this enhanced rate, notice 
that naively combining the data in Table 1-3 one obtains ^(ggF + ttH, 77) = 1.05 it 0.28 
and )u(VBF + VH, 77) = 1.8 it 0.6. It turns out that these different values are then realized 
with a slighly reduced -RggF and an increased R^^. 

The PDF of (3^ is asymmetric with a longer right tail. (3f, appears mainly in the h ^ bb 
decay rate, i.e. in Rj^^. ^^ The reason of the asymmetry is the following: as the branching 
ratio B{h — >• bb) = 57% in the SM, a deviation of /3fc from (hence Cj, from 1) results in a 
sizeable modification of the total width of the Higgs. Our signal strengths are expressed 
as jl{X,bb) = RxR^i/ Rmdth and contain a strong correlation between R^^i and i?width- As 



^h quark contributions in the loop-induced processes (ggF, 77, Z7) are small and can be disregarded. 
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Figure 5: On the left, posterior PDF of R^^ in scenario I (black) and scenario II (red). Also shown 
are the 2D posterior PDFs of i?width versus R-^-y (middle) and i?ggF versus i?^^ (right) in scenario 
I. Color code as in the previous figure. 




Figure 6: On the left, posterior PDF of i?width — r^/F™ in scenario I (black) and scenario II 
(red) . Also shown is the 2D posterior PDF of i?width versus i?ggF in scenario I (middle) and scenario 
II (right). Color code as in the previous figures. 

-Rwidth significantly increases with iJ^g, the deviations from //^^ = 1 are smaller than what 
we could naively expect, allowing large values of i?^^, hence /3b. This explains the tails of 
the PDF of /3b. 

The ID and 2D PDFs of i?width ^-re shown in the left pannel of Fig. 6. It turns out that 
a large increase of -Rwidth is not forbidden by the measurements of other channels, in which 
this effect is compensated by an increase of the decay or production rates, in particular 
ggF. The upper bound on i?width, Rmdth ^ 3, comes from the requirement /3{, < 1. 

In Fig. 7, we show the PDFs of the tensorial couplings Q.^, Qz-y^ Qz in Fig- 7 for scenario 
I. The PDF of C7 is constrained to small values in order to have the correct -ff — )• 77 rate. 
The PDF of Z7 is much broader because of the weak experimental sensitivity to the Z7 
rate. The distribution for Qz (and similarly Cv^) is mainly due to indirect effects on the 
fundamental parameters Pyy (77 and Z^ rates, as well as TGVs) rather than because 
of direct experimental constraints. Notice that even with the assumption of custodial 
symmetry (which enforces aw = az), Eq. (5.4) allows the rates for associated production 
to be different for Z and W because of the contribution of the tensorial couplings. It turns 
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Figure 7: Posterior PDFs of C^, ^27, and C^ in scenario I (black) and scenario II (red). 
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Figure 8: Posterior PDF of Rzu versus i?wH in scenario I (left) and scenario II (right). Color 
code as in the previous figures. 

out that C,w and C,z can be large enough in scenario I to induce a substantial deviation 
from i^wH = -RzH- This is shown in Fig. 8. This effect is also present in scenario II to a 
lesser extent. 

In scenario I, we observe two peaks of opposite signs for the tensorial couplings C,g 
and Q-f. These features appear because of the competition between the tree-level Cg,7 
and the loop-level SM couplings in the ggF and digamma amplitudes. In addition to the 
classical region where C, adds up to the SM coupling and cannot be very large, regions with 
C = 0(— 2A^^) are also allowed. Note that C7 is a linear combination of Pww, PwB, Pbb, 
but that the two ^^ peaks cannot be seen in the PDFs of these parameters. These four 
regions do not show up in scenario II, because the Q are loop-suppressed and thus cannot 
be large enough to cancel the SM couplings. 

A feature of the PDF of (zj is that, in spite of the various constraints on Pww, f^WB, 
and f^BB, the Z^ rate is still likely to be considerably enhanced. The shape of the (z-y PDF 
is mainly constrained by the CMS bound flz-y < 9.3 in scenario I, while indirect constraints 
from the S parameter and trilinear gauge vertices dominate in scenario II. In the enhanced 
rate, the HDO contribution dominates, such that Rzy is mostly proportional to {Cz-f)'^- 
This happens in scenario I, but also in scenario II although the Q are smaller. The PDF 
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Figure 9: Posterior PDFs of Rz^ = T{h -^ Z-f)/T'^^{h -^ Z-f) (left) and fiz^y (right) in scenario I 
(black) and scenario II (red). 

of the ratio Rz-y can be seen in Fig. 9. The 95% Bayesian credible intervals are [0, 12.0] for 
scenario I and [0, 4.3] for scenario II. As large deviations are allowed in this channel within 
this framework, it is therefore particularly promising for the discovery of a NP signal. 

In scenario I, Cz'y is sufficiently large to cancel the SM coupling, such that enhancement 
with both signs of Cz-y is realized. In contrast, for scenario II, only the branch with 
constructive interference (z'y < can enhance Rz'y- In both scenarios, (z-y can cancel the 
SM coupling such that having a small or vanishing Rz-f is likely. ^^ 

Finally, we compute the signal strength of /i — ?■ Z'j in case of a fully inclusive analysis 
at the LHC. The PDFs are shown in the right panel of Fig. 9 for both scenarios. In scenario 
I, the distribution reaches the CMS 95% C.L. bound jlz-y < 9.3, while it vanishes before in 
scenario II. The 68% and 95% BCIs are [0,3.6], [0,8.1] in scenario I and [0, 1.6], [0,3.2] in 
scenario II. 

Given that the 13/14 TeV LHC has a good potential to constrain the h — ?■ Z'j rate, one 
may wonder about the impact of a more precise measurement on our results. Therefore, 
we investigate the possibility of having flz-y < 2 at 95% CL, and we implement this bound 
as a step function. ^"^ It mainly results in a better determination of I3bb and /3\yw in 
both scenarios, as can be seen in Fig. 10. This new limit has an effect on the Higgs 
phenomenology in scenario I only. It leads to a slightly better prediction of -Rwh: the 
95% BCI is [0.7,1.5], instead of [0.5,1.7] when we take into account the current limit on 
h^ Zj. 



^^ We do not focus on this aspect as the direct searches at the LHC are still far from this level of precision. 
The shape of the Rzj PDF follows the distribution oi\(^z-y+^z-f\'^, which presents a peak in 0. Schematically, 

— 1/2 

for a uniform distribution of C,z-y, the peak behaves as ("^ . One can observe a similar behaviour for -Rttn- 
^^Note that the relative SM production rates c^x/^x '^x do not change significantly between 8 TeV 
and 14 TeV. Thus, assuming a fully inclusive analysis, we can take the decomposition into production 
channels as given in Sec. 3, Table 2. 
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Figure 10: Posterior PDF of f3ww versus /3bb in scenario I (left) and scenario II (right). The red 
and blue regions correspond to the 68% and 95% BCRs from the current measurements, while the 
black and grey contours correspond to the 68% and 95% BCRs assuming in addition that fiz-y < 2. 



8 Conclusion 

Less than one year after the first announcement of a signal in LHC data, the existence 
of a Higgs boson is at present firmly established. Time has come to probe in detail the 
structure of the now complete electroweak sector, searching for indirect signs of high-energy 
New Physics. 

We use a complete basis of dimension-six operators encoding NP effects in an effective 
Lagrangian in which all tensorial couplings are taken into account. The basis is chosen such 
that field-strength-Higgs operators (Off) are exactly mapped into tensorial couplings. In 
this basis it is straightforward to study the well-motivated hypothesis of loop-suppression 
of these operators. 

The data taken into account in our analysis are the whole set of results from ATLAS 
and CMS, including all available correlations, as well as Tevatron data. Trilinear gauge ver- 
tices measurements and constraints on electroweak precision observables are also included 
in our study. 

It turns out that weak bosons tensorial couplings can in principle modify non-trivially 
the kinematic structure of the VBF, VH and h — )• VV processes, and thus the efficiency 
of the kinematic selections. Scrutinizing the experimental analyses, we find that one can 
consider unchanged kinematic cuts efficiencies for the rates of these processes to a good 
approximation. The effect of tensorial couplings is crucial in the ggF, h — )■ 77, h — t- Z7 
processes, where they compete with the one-loop SM couplings. In our predicted rates, 
leading loop corrections are taken into account, consistently with respect to the loop- level 
HDO framework. 

In order to put constraints on the higher dimensional operators we consider, we carry 
out a global analysis in the framework of Bayesian inference. We find this approach par- 
ticularly appropriate as it deals with weakly constrained problems and naturally takes into 
account fine-tuning, such that the results we present are free of improbable, i.e. fine-tuned, 
cancellations. Markov Chain Monte Carlo techniques are used to perform the numerical 



28 



integrations. 

Our analysis is centered on New Physics arising at < 3 TeV. It allows us to express our 
results in terms of /3j = Oj v^/A^, i.e. to factorize the effect of the coefficients a, and of the 
scale of New Physics A. Among other things, this parameterization is A-independent up 
to a small log A dependence. Distributions of HDO coefficients at any low A can be easily 
mapped from our results. Two general scenarios are considered: I) democratic HDOs, 
where the coefficients of all the dimension-six operators are essentially unconstrained, and 
II) loop-suppressed Oj^f's, where the field-strength-Higgs operators are loop suppressed 
with respect to the other HDOs. 

We find overall a substantial amount of freedom in both of these scenarios. For in- 
stance, the coupling to the top quark still allows for 0{1) deviations, while the couplings 
to bottom and tau are slightly more constrained. We report in both scenarios a large and 
natural region at small or vanishing ct, favored by goodness of fit, which we trace back to 
the slight deficit observed in the bb channel with the ttH production mode, as reported 
by CMS. Regarding the couplings of the vector bosons, we find that only small deviations 
are allowed, of the order of 10% — 20%. Despite the fact that the inclusive rate of the 
/i — ;■ 77 channel is close to one, we observe a slightly enhanced Ryy = T^^/T?^^ , with the 
95% Bayesian credible interval (BCI) for R^y given by [0.8,2.5] ([0.8,2.3]) in scenario I 
(II). Correlations with -Rwidth = ^h/^f^ as well as the deviations of the various production 
modes from the SM are then identified to be responsible for fitting the experimental signal 
strengths. We also find that an enhanced coupling to the b quark is likely, leading to a larger 
total width, with the 95% BCI for i?width given by [0.7,2.7] ([0.6,2.5]) in scenario I (II). 
A strong correlation with RggF exists and makes the predicted signal strengths compatible 
with the data. 

Overall, it appears that the tensorial couplings play a crucial role in the SM loop- 
induced processes. In particular, after taking into account all the constraints we find that 
the Higgs boson decay width into Z7 can be enhanced by up to a factor 12 (4) in scenario 
I (II), within 95% BCI. Conversely, future measurements in the Zj channel will provide 
important bounds on the coefficients of the higher dimensional operators we consider. 
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Appendix 

A Deriving the tree-level effective Lagrangian 

Here we give some details on the derivation of the effective Lagrangian L^j + Ct- 
We note that the following field rescalings need to be performed in the SM Lagrangian in 
order to pass to canonically normalized kinetic terms for all fields: 



h 


-)• 


[l+Ch]/i, 




G, 


— ^ 


[1+C,]G^, 




K 


— )- 


[1 + i^\ A^ + 


iz-fZii , 


Z, 


— )- 


\\^iz\Z^, 






^ 


\^^iw\w^ 


1 



(A.l) 



with 



|«D2 + -a^2 ) V 



acG v^ 



?7 = I avKVK-s^ + olbbCw — -^ol-wbSwCw I 



t}2 



CZ7 = ( lOLWWCujSw — '^ObBCwSw — -^^awBic^u — S. ^ '' 
iz = f OLWWcij + OLBBsi, + l^aWBSwCw \ v^ , 

Cw = awwv^ ■ (A. 2) 

The quantities g, g' , v are replaced by the quantities g, Sw, v, which by definition are 
related to the fine-structure constant a, the Z mass mz and the Fermi constant Gp via 

|;2 2 -[^ 

47ra = slg^, m| = — ^ , Gf = ^=— • (A.3) 

Taking into account the HDOs we can compute in terms of the SM parameters 

ina = s^5^(l + 2sl,aww v^ + IcIjUbb v^ - CwS^awB v^) , 
2 v'^9'^ /^i , 1 ~2 , 1 / ~2 



z = -T=2~ [ 1 + 2"^^^ + 2"^^^ ~^ 2c^awwv + 2s^aBBV + CwSwawBv' 
Gf = -^(l-^an2V^ + ^aDvA . (A.4) 

For a, all corrections come from wave function renormalization (WFR) of Af^. For mz, 
there are both WFR and direct mass corrections, while for Gp = 5^/(4\/27tt-^) the WFR 



*The derivation follows closely Ref. [22]. 
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for W drops out and we are left with mass corrections as well as explicit vertex corrections. 
Combining Eq. (A. 3) and (A. 4) and inverting we obtain 

tP' = v^ - — (a£,2 - an) , 

g = g +g V - _ [a^2 + odJ " n _ n ^wB - 2aww , 

~2 2|222( 1 r'l li ^^wSw „r i\ 

^w = Sw + SujCujV [ „. 2 _ 2^ P-D2 + 0!d\ + ^ _ 2 "^^ + ^[OiWW " "BbJ I ■ 

(A.5) 

Next we need to express the bare Yukawa couplings in terms of the physical masses. One 

finds: 

\f2mf ( \ 2 -iX 

Vf = ( 1 + T i"D2 -anlv +afv 1 . (A.6) 

Finally we define a gs by A-kUs = g1 and write 

~gl = gl{l-2aGGv''). (A.7) 

B Loop-induced SM couplings 

Here we give the expressions of the form factors for SM loop- induced couplings Eqs. (2.19), 
(2.20). 

^/(r) = ^r(l + (l-r)/(r)), (B.l) 

^,(r) = ^(2 + 3T + 3r(2-r)/(r)), (B.2) 

Az,{t, k) = -^^(f{r) - f{K)) , (B.3) 

[T — K) 



(r — K) [T — ny (r — k,)^ 



Bz,{r, k) = j^—^ + -^—-^^{f{r) - /(k)) + j-—^M^) " 5('^)) • (B.4) 



arcsm^fr ^/^) if r > 1 

4( 



'*^'^^ i(.og|±^-..)Mf.;i. («-^' 



f \/7^^arcsin(r-i/2) if r > 1 , 
( 2 V 1 - T(log ^_^^_^ - ivr) if r < 1 . 

C Tensorial Higgs decay to weak bosons 

In the vectorial h — )• VV decay we encounter the following integral [64] 

Rs{x)^ j ^ dy^^— ^f-(y2-x) + l + x-2yj , (C.l) 
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where y = Ey/mh is the energy of the on-sheh gauge boson, x = rriy/mf^ and 6 = 0(1) for 
transversely (longitudinally) polarized gauge boson. In terms of the function 

F{x) = -=L= arccos (^-^) (C.2) 

^/4:X-l \ X2 J 

we can evaluate 

Ri + 2Ro = ^ f (-1 + x)(2 - 13x + 47a;2) 

— 3x(l — 6x + 4x ) log(x) 

+6x{l -8x + 20x)F{x) ) . (C.3) 



This is indeed proportional to the total decay width given in Ref. [56]. The translation of 
notation is 

Rt = l6{Ri + 2Ro) , (C.4) 

with the decay width given by 

r{h ^ VV*) = -^l-^mhdvRT{x) , (C.5) 

and Sw = I, Sz = j2- f^w + f^w- 

To compute the tensorial contributions, we need to multiply the integrand with uw = 
(y — x) and u? 



VV 



x + l 
2 



Q5ix)= I ~ dy ^^ _ ^Jl (^{y^-x) + l + x-2y]{y-x), 



x + l 



1 



Ps{x) ^ J ^ dy ^Y3^ (^- (/ - x) + 1 + X - 2yj {y - x)^ . 
The expressions are of similar form but slightly more complicated 

Qi + 2Qo = -^ U-1 + x)(-3 + 31x - 140x2 + leOx^) 

+ 3x(2 - 15x + 42x2 - 12x^) log(x) 

+6x(-2 + 19x - 68x2 _^ s4:X^)F{x) J , (C.6) 

p^ + 2Po = ^^^ ( (-1 + a;)(12 - 173x + 971x2 _ 2299x3 + 1405x^) 
- 6x(5 - 46x + 174x2 - 288x3 + 48x*^) log(x) 

+12x(5 - 56x + 256x2 - 544x3 + 432x'')F(x) J . (C.7) 

The final quantities v'vVj ^vy introduced in Eq. (4.7) are 

Qi + 2Qo / 2 \ Pi + 2Po /^ gx 
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